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The Garvey-Kelson (GK) relations are powerful algebraic expressions connecting the masses of 
neighboring atomic nuclei and derived under reasonable physical assumptions. In this contribution 
we show that these relations are even more general than originally assumed as their validity is 
conditioned by the smoothness of the underlying nuclear mass function and nothing else. Based 
on the assumed smoothness of the underlying physical function, the main conclusions from this 
study are: (1) the GK relations are model independent; (2) any slowly- varying physical observable 
satisfies the GK relations; and (3) the accuracy of the GK relations may be systematically improved. 
Examples from nuclear physics (charge radii) and from particle physics (baryon masses) are used to 
illustrate the flexibility of the approach. 
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Nuclear observables display a remarkably smooth de- 
pendence on atomic number (Z), neutron number (TV), 
and baryon number (A = Z+N). Perhaps the best known 
example of such a smooth behavior is the nuclear mass 
formula of Weizsacker [1 and Bethe and Bacher [2]. In 
a formulation that dates back to 1935-36 and that has 
survived the test of time, the atomic nucleus is mod- 
elled as an incompressible liquid drop consisting of two 
strongly interacting quantum fluids (one neutral and one 
charged). The mass formula includes volume, surface, 
Coulomb, and asymmetry terms that vary smoothly with 
A and Z. When these four empirical coefficients are fit- 
ted to the enormous database of ground-state masses (of 
more than 2000 nuclei [3 ) a rms deviation of about 3 
MeV is obtained. Most of the 3 MeV discrepancy can be 
explained by invoking physics that ceases to be smooth, 
such as pairing correlations and shell corrections. Indeed, 
highly successful modern mass formulas that incorporate 
these effects are able to reduce the rms deviation by al- 
most an order of magnitude, to about 300 keV [4j|5]. 

An approach that has been recently revitalized, both 
because of an interest in understanding any inherent lim- 
itation in the nuclear-mass models as well as because of 
its astrophysical applications, is the one by Garvey and 
Kelson E [TO] . The Garvey-Kelson (GK) relations 
are powerful algebraic expressions connecting masses of 
neighboring nuclei. Based on a few simple assumptions 
about the underlying nuclear dynamics, Garvey, Kelson, 
and collaborators derived the following two relations con- 
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necting the masses M of six neighboring nuclei [6j [7] : 

AM 6 (1) (TV, Z) = M(N + 2, Z - 2) - M(N, Z) 

+ M(N, Z - 1) - M(N + 1, Z - 2) 

+ M(N + 1, Z) - M(N + 2, Z - 1) = , 

(la) 

AM 6 (2) (TV, Z) = M(N + 2, Z) - M(N, Z - 2) 

+ M(N + 1, Z - 2) - M(N + 2,Z-1) 
+ M(N, Z - 1) - M(N + 1, Z) = . 

(lb) 

It is pertinent to note that in their original publication 
Garvey and Kelson made the following germane state- 
ment: "No assumptions are made about the quantitative 
aspects of the description except that the position of the 
single-particles levels, and the residual interactions be- 
tween nucleons in them, vary slowly with atomic num- 
ber" [6]. In some sense, this statement is reminiscent of 
the role that symmetries and group theory have played in 
Nuclear Physics and that had the late Marcos Moshin- 
sky as one of its leading exponents (see, for example, 
Ref. [11). Rather than attempting an exact solution of 
the difficult many-body problem, group-theoretical ap- 
proaches assume an underlying symmetry of the (often 
unknown) Hamiltonian. Based on such an assumed sym- 
metry, powerful relations can then be derived among 
the various states of the Hamiltonian that are connected 
through the symmetry. It is the aim of the present Letter 
to explore the consequences of the slowly-varying dynam- 
ics proposed by Garvey and Kelson and to establish the 
power and generality of the GK relations, assuming only 
the smoothness of the underlying physical observables. 

In particular, it will be shown that the linear combina- 
tions of masses chosen by Garvey and Kelson are propor- 
tional to the third-order derivative of the underlying nu- 
clear mass function M(N,Z) |l2 r That is, the selected 
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combinations of masses are insensitive to the underlying 
mass function as well as to its first and second deriva- 
tives. It is evident then, that as long as the mass func- 
tion is slowly varying, the GK relations will be satisfied 
to a very good approximation. For example, if M(7V, Z) 
is represented by the smooth liquid-drop formula, then 
the GK mass relations are satisfied to better than 150 
keV for A > 50 and to better than 50 keV for A > 100, 
for nuclei close to the line of stability. When using the 
available experimental database of nuclear masses, Barea 
and collaborators [HI [H [10] have shown that the GK re- 
lations are accurately satisfied throughout the periodic 
table (on average, to about 180-200 keV). Thus, it was 
concluded that the GK relations may provide a unique 
tool to test and improve current nuclear mass formu- 
las [H [10] — especially in the as yet experimentally un- 
explored regions of exotic nuclei. 

To express the GK relations in terms of the third-order 
derivatives of the nuclear mass function M(N,Z), both 
N and Z are treated as continuous variables so that a 
Taylor series expansion around the reference point (TV, Z) 
may be performed. It is then a simple algebraic exercise 
to show that the two GK relations are identically equal 
to: 

<^™- + < 2b » 

Whereas Garvey and Kelson invoked reasonable physi- 
cal arguments to construct the two linear combinations 
of masses given in Eq. ([!]) , the previous result shows that 
the GK relations are both more general and more pow- 
erful than originally assumed. Indeed, the GK relations 
are model independent as they are insensitive to the un- 
derlying dynamics, provided such dynamics generates a 
slowly- varying mass function M(N,Z). Moreover, we 
observe that any smooth physical observable — not only 
M(N, Z)— will satisfy the GK relations. 

As noted in Eq. (J2|, the GK relations are proportional 
to the crossed (or mixed) derivatives of M(N,Z). As 
any function of two variables has four independent third- 
order derivatives, two more "Garvey-Kelson" relations 
may be derived to this order. That is, 

AM^(N, Z) = M(N + 2,Z)- 3M(N + 1, Z) 



•3M(iV, Z) - M(N - 1, Z) 
d 3 M 



dN 3 



0(<9 4 M) , 



AM| 2) (iV, Z) = M(N, Z + 2)- 3M(iV, Z+l) 
+ 3M(N, Z) - M(N, Z-l) 
d 3 M 



dZ 3 



0(<9 4 M) . 



(3a) 



(3b) 



N) chain, respectively. Note, however, that the formal- 
ism is completely general so one is free to select any two 
independent variables (e.g., A and N — Z) to carry out 
the analysis. Finally, note that the approach may be 
systematically improved by selecting linear combinations 
of masses so that third-order (and even higher-order) 
derivatives get cancelled out. Such extended GK rela- 
tions afford higher accuracy in the prediction of nuclear 
masses, albeit at the expense of having reliable infor- 
mation on the masses of a larger number of neighboring 
nuclei. Indeed, Barea and collaborators have recently ob- 
tained a single relation involving the masses of 21 nuclei; 
see Eq. (4) of Ref . [9] . If one denotes such linear combi- 
nation as AM2i(Af, Z), one can show that all derivatives 
up to fifth order get cancelled out! That is, 



AM 21 (iV, Z) = 2 



d 6 M 



d 6 M 



dZ 2 dN A dZ A dN 2 



■0(d 7 M) . (4) 



Thus, one can systematically improve the accuracy of the 
original GK relations by progressively removing third- 
and higher-order derivatives. 

We next present two examples that aim to illustrate 
the generality and flexibility of the approach. One of 
them goes beyond nuclear masses and tests the validity 
of the GK relations for nuclear radii. The other one is 
used to derive — or rather to re-derive — some well-known 
mass relations among baryons. The chosen examples are 
not intended to be exhaustive or systematic, but serve to 
illustrate the potential of the GK relations. 

Our first example concerns itself with nuclear radii. 
The charge distribution of atomic nuclei can be accessed 
very cleanly with electromagnetic probes and, indeed, the 
charge radii of many nuclei are known to exquisite ac- 
curacy. Although smaller than for masses, there exists 
a sizable experimental database of nuclear charge radii 
(see the 2004 compilation by Angeli [33]). In applying 
the GK relations to nuclear charge radii, we find a rms 
deviation of about 0.01 fm for the approximately 130 nu- 
clei that satisfy (at least) one of the two original GK 



TABLE I: Central value of the measured charge radii for 7 
of 8 nuclei required to test the two Garvey-Kelson relations 
(l) in the region around 195 Pb (Z = 82, iV=113), taken from 
Ref. [13] . The charge radius of 195 Pb is unknown. Theoretical 
neutron radii obtained with the FSUGold interaction [14] are 
also shown. 



These last two relations correlate the masses of four nu- 
clei along an isotopic (constant Z) and isotonic (constant 
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19 slPb 


5.4420 


5.5390 


^Hg 


5.4347 


5.5355 




5.4239 


5.5102 






5.5143 


80 ttg 


5.4311 


5.5219 


194IT11 
81 11 


5.4233 


5.5119 


196rpi 

81 11 


5.4304 


5.5369 


^Hg 


5.4311 


5.5219 


194rpi 

81 11 


5.4233 


5.5119 




5.4420 


5.5259 


x gpb 


5.4420 


5.5259 


196rpi 

81 11 


5.4304 


5.5369 


\AR^\ 


0.0001 


0.0002 


\AR^\ 


5.4385 


0.0002 



3 



relations [see Eq. Q]. This result is consistent with the 
smooth global structure of the charge radius identified 
in Ref. [13]. To provide a specific example, we list in 
Table [T] the charge radii of 7 out of 8 nuclei that are re- 
quired to test the GK relations ([!]) in the 195 Pb region. 
A pictorial representation is also provided in Fig. [I] The 
available charge radii are sufficient to verify that the GK 
relation ( lb ) is satisfied — at least in this region — to a re- 
markably high accuracy (~10 -4 fm). Having established 
the validity of one of the GK relations in this region, 
one can then use the other relation to predict the central 
value of the charge radius of 195 Pb which is unknown ex- 
perimentally to be i? ch ( 195 Pb) = 5.4385 fm. Note, how- 
ever, that the GK relations ([!]) applied to 195 Pb with 
different neighboring nuclei allow up to 12 independent 
estimates [9] for i? c h( 195 Pb). For the present case, the 
experimentally available data limit these to four. When 
all these four estimates are averaged, one obtains a value 
of # c h( 195 Pb) = 5.437(5) fm. For completeness, theoret- 
ically generated neutron radii in the same region have 
been added to the table. Although unavailable experi- 
mentally, the close fulfilment of the GK relations suggests 
that the underlying neutron-radius function R n (TV, Z) 
(whatever it may be) is likely to be smooth. 

In our second example we derive Garvey-Kelson re- 
lations for baryon masses. For simplicity, we focus ex- 
clusively on the baryon octet and decuplet. To apply 
the GK relation to baryons, it is sufficient to replace the 
neutron and proton numbers in Eqs. ([I]) and ^ with the 
strangeness (N — > S) and electric charge (Z —> Q) of the 
baryon, respectively, and to expand around S = — 2 and 
Q = 1. As there is no double-strange baryon of positive 



charge, the GK relation (la) is not applicable. However, 
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FIG. 1: Pictorial representation of the Garvey-Kelson rela- 
tions in the 195 Pb region [Z = 82, N = 113). The unshaded 
(shaded) circles represent six nuclei involved in the first (sec- 
ond) Garvey-Kelson relation. The "plus" and "minus" signs 
inside the circles indicate how they should be included in the 
GK relations. The nucleus 195 Tl plays no role in this example. 



the relation ( lb ) may be applied and yields the interest- 
ing result 



{p-n) 
d 3 mg 



d 3 ms 



dQ 2 dS dQdS 2 



) + (s°-s-) 

0(d 4 m 8 ) , 



(5) 



where mg(S,Q) is the underlying ground-state baryon 
octet mass function. The equation Acg = is the cel- 
ebrated Coleman- Glashow (CG) mass relation [15 ; de- 
rived originally using unbroken flavor SU(3)[15] and later 
also by methods such as the 1/N C expansion [16j[T7J[T8] 
and the QCD parametrization method [19] . The CG mas 
relation vanishes within experimental error [I6j [17| [THJ 
[19] • This suggests that the underlying mass function 
m$(S, Q) — although likely unknown — is smooth enough 
to justify the use of the GK relations in the baryon sec- 
tor. However, to this order, the GK relations yield no 
further information because the structure of the octet is 
too limited for the application of Eq. |3| . 

Although the first of the GK relations is still not appli- 
cable in the case of the (J 71 = 3/2 + ) baryon decuplet, the 
three other relations ( |lb| ), ( |3a| ), and ( |3b| ) yield interesting 
results, namely, 

A£ G = (A+ - A ) - (£*+ - £*-) + (~*° - £*-) * , 

(6a) 

A Q= _i = A" - 3S*- + 3~*~ - O « , (6b) 
A s= = A++ - 3A+ + 3A° - A" « . (6c) 

Note that the approximate equal sign implies that all 
three expressions are proportional to a third derivative of 
an underlying (decuplet) mass function that is assumed 
to be slowly varying. Indeed, all three equations vanish 
to within experimental errors. The first of the three equa- 
tions is the corresponding Coleman-Glashow relation ex- 
tended to the baryon decuplet. The second u iso charge" 
(Q = — 1) relation connects the masses of all negatively- 
charged baryons in the decuplet. Finally, the last "isos- 
trangeness" (S = 0) relation connects the masses of the 
four A isobars in the multiplet (this relation is referred as 
A3 in Ref. [E]). Note that there is an additional relation 
connecting four of the baryons in the decuplet, vis-a-vis, 



±Q-S=2 



A++ - 3S 



*+ 



3E 







(?) 



however, it is not independent of the other three. 

It should be noted that essential to the derivation of 
the mass relations is only the structure of the multiplets, 
namely octet and decuplet, but not the particular use 
of the ground-state baryons. Indeed, identical mass re- 
lations will hold true whether the hadrons are baryons 
or mesons or whether they contain strange or charmed 
quarks. Thus, we conclude this section with a remark 
on the masses of charmed baryons. Replacing all strange 
quarks with charm quarks in Eq. ^ yields the following 
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(J 71 " = 3/2 + ) mass relation: 

AgL 2 = A++ - 3E++ + 3H++ - n++ « . (8) 

Since at present neither the mass of the (3/2 + ) double- 
and triple-charm baryons are known, a prediction can 
be made only for a specific linear combination. Using 
the most up to date information from the Particle Data 
Group [20] we obtain 

3H++ - ntct ~ - A++ « 6328 MeV . (9) 

Note that a recent theoretical study based on the Bethe- 
Salpeter equation yields a value of 6396 MeV for this 
combination [21 j, or a 1% discrepancy. 

In summary, although well motivated on physical 
grounds, we found the Garvey-Kelson relations signifi- 
cantly more general and more powerful than originally 
assumed. Using the simple — yet little-known — fact that 
the linear combinations of masses that enter into the GK 
relations are proportional to the third derivatives of the 
nuclear mass function [12] , we concluded that the validity 
of the GK relations hinges exclusively on the smoothness 
of the underlying function and on nothing else. Having 
established this fact, we regard the following three as the 
main conclusions of this work: (1) the GK relations are 
model independent; (2) any slowly-varying physical ob- 
servable satisfies the GK relations; and (3) the accuracy 
of the GK relations may be systematically improved. 

As examples chosen to illustrate the power and flexi- 
bility of the approach, the GK relations were extended to 
the calculation of nuclear radii and mass relations in the 
baryon sector. In the case of nuclear radii, we focused on 
a region of the periodic table that enabled both the vali- 
dation of the approach and the prediction of the unknown 
charge radius of 195 Pb. For the case of baryon masses, 
a direct application of the GK relations to the ground- 
state octet and decuplet reproduced some well-known 
mass relations, such as the celebrated one by Coleman 
and Glashow [15 . In addition, a model-independent pre- 
diction was made for a linear combination of the masses 
of the double- and triple-charm baryons that seems to be 
in good agreement with available theoretical calculations. 

In the future we plan to a carry out a comprehensive 
and systematic program that will extend the Garvey- 
Kelson approach to areas that go beyond its original in- 
tent. However, we hope that the present study may al- 
ready inspire the reader to consider applications in other 
fields. We are confident that the approach can be success- 
fully generalized to the study of physical observables in 
atomic, nuclear, and particle physics — and perhaps even 
to areas outside of physics. The sole requirement for 
the Garvey-Kelson relations to hold is a slowly-varying 
observable. Remarkably, nature often seems to gener- 
ate precisely such kind of observables independent of the 
complicated — and likely unknown — underlying dynam- 
ics. 
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